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ABSTRACT 

Results  are  presented  from  the  study  of  the  off-axis  X 
and  Gamma  radiation  field  caused  by  a  highly  relativistic 
electron  team  in  liquid  Nitrogen  at  various  path  lengths  out 
to  2  radiation  lengths.  The  cff-axis  dose  in  Silicon  was 
calculated  using  the  electron/photon  transport  code  CY1TRAN 
and  measured  using  thermal  luminescent  dosimeters  (TLD' s)  . 
Calculations  were  performed  on  a  CdC-7600  computer  at  Los 
Alamos  National  Laboratory  and  measurements  were  made  using 
the  Naval  Postgraduate  School  100  Mev  Linac.  Comparison  of 
the  results  is  made  and  CYLTRAN  is  found  to  be  in  agreement 
with  experimentally  measured  values.  The  CYLTRAN  results  are 
extended  to  the  off-axis  dose  caused  by  a  100  Mev  electron 
beam    in    air   at    Standard    Tempature   and   Pressure    (STP)  . 
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I.     INTRODUCTION 

A.       BACKGROUND 

When  a  high  energy  electron  passes  through  matter  it 
loses  energy  through  two  classes  of  interactions;  ionization 
and  bremsstrahlung.  At  lower  energy  ionization  is  the 
predominant  energy  loss  mechanism  while  at  higher  energy 
bremsstrahlung  becomes  predominant.  Equation  1.1  shows  the 
energy  loss  by  a  beam  of  highly  relativistic  electrons  due 
to  ionization.  Equation  1.2  shows  the  energy  loss  due  to 
bremsstrahlung    for    a   highly    relativistic    electron.      [Ref.    1] 

(d2/dx)ion   =  (UPie4  n/mc2  )  (log  (  2mc2  /I)  (eqn    1.1) 

-    3/2(log(1-B    ))     -    1/2(log3)    +    1/16} 

<dE/dx>rad=    (4NZ2  R2,  hk/1  37)   [log  (  183/Z,/3)  }  (eqn    1.2) 

The  approximate  ratio  of  energy  loss  as  a  function  of 
distance  for  the  two  mechanisms  is  shown  in  equation  1.3 
where  E  is  the  electron  energy  in  JleV.  Table  I  contains  the 
meaning    of    the    symbols    used    in    these    equations 

In  the  case  where  the  product,  EZ  is  much  greater  than 
800  MeV,  energy  loss  by  bremsstrahlung  predominates  and  the 
energy  loss  due  to  ionization  may  be  neglected.  The  energy 
at  which  the  ratio  is  equal  to  one,  for  a  particular  Z ,  is 
known  as  the  critical  energy  and  is  denoted  E (crit) .  The 
critical  energy  in  Mev  is  shown  in  equation  1.4  .  If,  in  a 
given  material,  the  incident  electron  energy  is  much  larger 
than  E(crit)  then  effectively  bremsstrahlung  is  the  only 
energy  loss  mechanise  until  the  energy  of  the  electron 
approaches    E  (crit) .       [Ref.     1] 


TABLE   I 
Symbols   used   in    Transport  Equations 

Symb 

2i 

Meaning 

Pi 

3.14.. . 

n 

Electron    density 

e 

Charge    on  electron 

m 

Ilass   of    electron 

c 

Speed    of    light 

I 

Ionization    potential 

E 

Ratio    of    electron    velocity    to    c 

h 

Planck's   constant 

N 

Atomic    density 

Z 

Atomic    number 

E 

Classical   radius   of   the   electron 

k 

„..  . 

Frequency   of    the    highest    energy    photon 

produced 

(dE/dx) rad/(dE/dx) icn=EZ/800  (egn    1.3) 

E (crit) =800/Z  (eqn    1.4) 

Photon  production  by  bremsst rahlung  is  essentially  a 
statistical  process  where  the  energy  of  any  given  photon 
produced  in  an  interaction  follows  a  probability  distribu- 
tion function  that  is  determined  by  the  energy  of  the  elec- 
tron and  the  atomic  number  of  the  material.  The  spectrum  of 
energies  produced  by  a  mono-energetic  electron  beam  will  be 
very  broad,  with  some  photons  having  a  maximum  energy  equal 
to    that    of    the    incident    electron. 


The  bremsstrahlung  produced  photons  undergo  two  types  of 
interactions,  Compton  scattering  and  Pair  Production.  At 
energies  above  about  10  MeV,  pair  production  dominates  the 
energy  loss.  The  electrons  produced  in  these  interactions 
and  the  positron  produced  in  pair  production  undergo  further 
bremsstrahlung .  The  pair  will  eventually  undergo  pair  anni- 
hilation and  produce  two  low  energy  gamma  rays.  This  process 
of  bremsstrahlung  followed  by  electron/positron  materializa- 
tion followed  by  further  bremsstrahlung  leads  to  the  forma- 
tion of  what  is  called  a  'radiation  shower'  [Ref.  2] 
[Ref.  3].  Figure  1.1  is  a  graphical  representation  of 
shower  development. 

The  work  previously  cited  fcy  Oppenhimer  and  later  work 
by  Rossi  was  intended  for  energies  greater  than  the  critical 
energy.  This  early  work  was  intended  to  explain  the  high 
energy  showers  detected  in  the  atmosphere  as  a  result  of 
cosmic  rays.  With  the  advent  of  electron  accelerators  in  the 
hundred  rteV  and  GeV  range,  their  work  was  applied  to  the 
problem  of  dose  buildup  in  shielding.  Since  most  shielding 
is  made  of  high  Z  material  the  approximation  that  ionization 
can  be  ignored  is  valid  [Ref.  4]. 

Early  work  on  radiation  showers  treated  energy  loss 
by  ionization  as  just  a  disappearance  of  energy  from  the 
beam,  and  did  not  attempt  to  treat  the  electrons  formed  by 
the  ionization.  Since  these  theories  ignore  the  ionization 
electrons  contributicn  to  the  total  number  of  electrons 
present  in  the  shower,  any  dose  calculation  based  on  these 
theories  would  be  in  error  by  an  amount  egual  to  the  contri- 
bution of  those  electrons.  For  showers  formed  by  extremely 
energetic  electrons  the  number  of  ionization  formed  elec- 
trons will  remain  negligible  for  a  considerable  distance 
[Ref.  3].  In  energy  regimes  where  the  incident  energy  is  on 
the  order  of  the  critical  energy,  about  102  MeV  for  air,  the 
error  in  a  computed  dose  can  be  very  large. 
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REPRODUCED    AT    GOVERNMENT    EXPENSE 


Figure  1.  1    Shower  Development. 
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In  practice,  it  becomes  necessary  to  use  Monte-Carlo 
techniques        to  calculate         the  shower-induced        dose. 

Historically,  these  codes  have  been  used  in  a  wide  variety 
of  problems  dealing  with  shielding  and  dose  buildup  [Ref.  5] 
[Ref.  6].  There  appears  little  published  work  comparing 
calculated   and    measured    dose    rates   in    low    Z    material. 

B.       DESCRIPTION 

It  is  the  purpose  of  this  thesis  to  determine  the  off- 
axis  radiation  dose  deposited  in  Silicon  by  a  100  MeV  elec- 
tron team  in  liquid  Nitrogen,  at  ranges  up  to  two  radiation 
lengths  and  extend  these  result  to  the  off-axis  dose  in  air 
at  SIP.  This  energy  is  nearly  e^ual  to  the  critical  energy 
which  is  about  102  MeV  [Ref.  4].  The  energy  of  an  electron 
beam  is  a  function  of  penetration  depth  and  the  atomic 
number  of  the  material-  In  the  case  of  a  mono-energetic 
electron    beam,    the    beam   energy    as   a    function   of    depth    is: 

<E>    =    EQ  exp-  (x/L)  (eqn    1.5) 

where  E  is  the  incident  energy,  x  is  the  depth  in  the 
material  and  L  is  the  radiation  length.  The  radiation  length 
is    [Ref.    1  ]    : 

1/L=(4Z2  NR^/137)  log  (133/Zl/3)  (eqn    1.6) 

The  radiation  length  in  the  standard  80%  nitrogen  20% 
oxygen  atmosphere  at  STP  (standard  temperature  and  pressure) 
is  about  307  meters.  This  poses  a  significant  problem  with 
the  Naval  Postgraduate  School  Electron  Linear  Accelerator 
(linac)  which  has  an  experimental  area  only  two  meters  in 
length.         The    radiation      length    for    liquid    Nitroqen      is    47. 7 
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centimeters.  Since  Nitrogen  makes  up  eighty  percent  of  the 
atmosphere  and  the  next  most  abundant  component,  Oxygen 
differs  by  an  atomic  number  of  one  from  Nitrogen,  the  inter- 
action cross  section  of  Nitrogen  closely  approximates  that 
of  the  atmosphere.  In  order  tc  measure  the  off  axis  radia- 
tion field  at  several  radiation  lengths  it  was  decided  to 
substitute  liquid  Nitrogen  for  air.  Calculations  where  made 
using  the  electron/photon  transport  code  CYLTRAN  and  are 
compared  to  measurements  at  26,52,78  and  104  centimeters  of 
liquid  nitrogen.  The  off  axis  radiation  field  was  calculated 
and  measured  at  1  centimeter  intervals  from  on-axis  to  ten 
centimeters  off-  axis.  The  lengths  used  in  both  the  calcu- 
lations and  experiments  were  26,52,78,104  centimeters  of 
liquid  nitrogen  corresponding  tc  0.54,  1.08,  1.63  and  2.17 
radiation  lengths  respectively,  Table  II  compares  the  exper- 
imental lengths  to  lengths  in  air  including  the  off-axis 
distances   covered. 
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The  assumption  that  off-axis  distances  scale  linearly 
with  radiation  length  is  based  on  the  work  of  Oppenhimer  and 
Rossi.  Their  work  shows  that  the  development  of  a  radiation 
shower  is  directly  related  to  the  position  measured  in  units 
of  the  radiation  length.  The  lateral  spread  of  the  shower 
is  also  dependant  on  the  distance  in  units  of  radiation 
length.  For  materials  with  similar  Z,  it  is  clear  from  equa- 
tion 1.6  that  the  radiation  lengths  are  in  proportion  tc  the 
ratio   of   the   atomic    densities   £  fief .    2]   [Ref.    3]. 
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II.  CALCULATIONS 

A.   DESCEIPTION 

The  cylindrical  geometry,  multi-material  electron/photon 
Monte-Carlo  transport  code  CYLTEAN  was  used  to  calculate  the 
off-axis  dose  deposited  in  a  0.2  cm  thick  silicon  target  at 
the  desired  ranges.  The  silicon  target  was  divided  into  ten 
concentric  rings  of  width  1  centimeter  centered  on  the  axis. 
The  simulation  geometry  is  shown  in  Figure  2.1  .  The  portion 
of  the  CYLTEAN  output  of  interest  to  this  experiment  is  the 
energy  deposition  in  each  material  zone.  This  portion  gives 
the  material,  its  location  and  mass,  and  the  energy  depos- 
ited by  primary  and  knock-on  electrons,  by  photons  and  the 
total  energy  deposited.  Each  value  for  the  energy  deposited 
includes  a  one  or  two  digit  number  that  represents  the  one- 
sigma  statistical  uncertainty  of  the  quantity.  The  deposi- 
tion data  is  produced  in  Mev  per  incident  electron.  To 
simplify  future  comparisons  this  data  was  converted  to  a 
normalized  dose  in  urits  of  F.ads  (Si)  per  coulomb  .•  Table 
III  contains  the  total  calculated  normalized  doses,  and 
Figure  2.2  through  Figure  2.5  are  graphs  of  the  off-axis 
results.  Figure  2.6  is  a  comparison  of  the  calculated  off- 
axis  normalized  doses.  Appendix  A  contains  the  complete 
CYLTEAN  energy  deposition  data  and  further  information  on 
CYLTEAN. 
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Figure  2.1   CYLTEAN  Problem  Geometry. 

17 


q; 

W 

o 

Q 

W 

•H 

H 

X 

H 

XI 

H 

1 

4-1 

pq 

4-1 

^ 

O 

pq 

<! 

T3 

Eh 

0) 

-P 

id 

rH 

3 

u 

M 

(0 

U 

— —  - —  — 

!» 

I"" 

— " 

*"" 

<~ 

vO 

1*- 

in 

m 

\o 

•o 

^o 

•" 

N 

s 

f*. 

N 

N 

N 

n 

N 

N 

N 

O 

o 

o 

O 

o 

o 

o 

O 

O 

o 

• 

T"" 

r— 

r~ 

r™ 

r- 

r- 

r— 

r— 

r- 

r- 

s 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

^-» 

u 

r^ 

rn 

on 

rn 

in 

in 

r^ 

=t 

co 

in 

5T 

VO 

rn 

CO 

^o 

VO 

o> 

CO 

o 

r« 

ja 

^r 

CO 

T— 

CO 

in 

-O 

m 

m 

co 

^o 

CN 

a 

o 

0 

t— 

LO 

r- 

m 

a- 

rn 

CN 

CN 

T— 

r- 

i— 

.H 

a 

o 

u 

rn 

o 

o 

o 

o 

CN 

r— 

T— 

r~ 

r- 

in 

r^ 

=* 

r- 

r- 

cn 

1 

C|U 

(Diaj 

CO 

CD 

CD 

N 

N 

N 

r- 

N 

N 

r- 

da. 

o 

o 

o 

O 

o 

o 

o 

o 

O 

O 

Ol 

1 

MU-. 

• 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

+J|H 

a 

r- 

«-" 

CO 

cn 

CN 

r-« 

cn 

-3" 

CO 

m 

•HtO 

u 

co 

CO 

en 

r- 

CO 

cn 

r- 

cf 

cn 

cn 

IzM 

CO 

r™ 

in 

r— 

co 

r- 

CO 

m 

r^ 

r- 

o 

T3lW 

r- 

cn 

r- 

t— 

co 

r» 

3- 

=j- 

m 

CN 

CN 

•Hta 

3I(T3 

ciw 

•Hi 

-lie 

r-» 

r^ 

r» 

CO 

in 

3- 

in 

r» 

m 

CO 

•H 

<+-ll 

CD 

CO 

CD 

co 

CO 

CO 

CO 

00 

K 

N 

Old) 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

W 

-ClO 

• 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

I     +jfa 

a 

r- 

rn 

J- 

o 

CN 

CT\ 

3- 

■=T 

r^ 

<J\ 

tTi 

0 

r^ 

rn 

CN 

00 

r- 

CN 

o 

O 

(N 

CN 

«ra 

r* 

U3 

CM 

rn 

r~ 

»"" 

r> 

CN 

r» 

rn 

(V|0J 

fN 

H-)lN 

ID 

vO 

LD 

^T 

m 

CN 

CN 

r™ 

<— 

cn 

r- 

•H 

H 

id 

s 

O 

vO 

m 

u 

=r 

ro 

r\! 

rn 

rn 

in 

CO 

r— 

r— 

<N 

0 

G 

0< 

o> 

o 

o> 

CD 

CO 

N 

N 

^ 

■0 

*""" 

O 

o 

o 

o 

o 

o 

O 

O 

o 

o 

• 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

a 

ro 

rn 

r* 

CM 

r» 

m 

in 

»— 

CO 

CO 

u 

r- 

r~ 

o 

o 

rn 

m 

m 

o 

v£> 

CO 

r» 

m 

T— 

=r 

(N 

rn 

CO 

r^ 

r» 

in 

^D 

CN 

r-« 

in 

n 

r~ 

in 

CN 

CO 

3" 

CN 

r- 

win 

o 

UIH 

o 

o 

o 

o 

o 

o 

o 

O 

o 

t 

C|X 

o 

rOlcO 

r— 

CN 

on 

C3- 

in 

U3 

["*■ 

00 

cn 

r- 

•PI 

1 

l 

1 

1 

1 

1 

1 

1 

1 

1 

WIH 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

■hw  a 

OIOU 

o 

fN 

m 

=t 

in 

lO 

r- 

CO 

cn 

REPRODUCED    AT    GOVERNMENT    EXPENSE 


0.0 


LEGEND 
E     CALCULATED  DOSE 


2.5 


5-0  7.5  10.0 

DISTANCE  OFF-AXIS  (CM) 


12.5 


Figure   2.2        Calculated   Off-axis   Normalized   Dose   at    26   ca. 
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CYLTEAN  is  a  relatively  large  simulation  code  of  about 
30,000  lines  of  fortran.  The  calculations  were  performed 
during  an  experience  tour  at  the  Los  Alamos  National 
Laboratories  on  a  CDC-7600  computer.  The  statistical  uncer- 
tainties can  be  seen  to  grow  as  the  thickness  of  the  liquid 
Nitrogen  is  increased-  This  is  due  to  the  exponential  rise 
in  the  number  of  interactions  that  the  program  must  simu- 
late. One  method  that  was  used  to  decrease  the  statistical 
uncertainties  was  to  increase  the  number  of  incident 
photons.  The  drawback  to  this  nethod  is  that  the  number  of 
interactions  has  already  increased  exponentially  as  the 
length  increases,  and  adding  additional  incident  electrons 
increases  the  cpu  time  of  the  program  even  further.  It  was 
found  that  an  increase  by  a  factor  of  two,  decreased  the 
statistical  uncertainties  by  only  ten  percent  but  increased 
the    cpu      time    by   at    least      a    factor    of    two.  At    twenty-six 

centimeters   and      with  five      thousand    incident      electrons    the 
cpu      time   was      five    minutes;  at    104      centimeters    and      with 

50,000    incident      electrons    the    cpu      time    increased      to    about 
two    hundred   minutes. 
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Figure    2.3        Calculated   Off-axis   Normalized   Dose   at    52   c; 
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Figure    2.4        Calculated   Off-ajcis    Normalized   Dose    at    78   ca. 
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Figure   2.5        Calculated   Off-axis   Normalized   Dose   at    104   cm. 
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III.  EXPERIMENTS 

A.   DISCUSSION 

All  experimental  runs  weie  conducted  at  the  Naval 
Postgraduate  School  using  the  Electron  Linear  Accelerator 
located  there.  A  full  description  of  tne  NPS  Linac  is  given 
in  [  Eef  .  7].  Two  parameters  required  close  monitoring 
during  experimental  runs;  the  electron  energy  (all  calcula- 
tions were  made  for  100  MeV)  and  the  total  charge  delivered 
by  the  accelerator. 

The  total  charge  delivered  vas  monitored  by  a  secondary 
emission  monitor  (SEM)  that  was  used  with  a  vibrating  reed 
electrometer  to  determine  the  tctal  charge  delivered  to  the 
SEM.  Once  the  total  charge  to  the  SEM  was  known,  the  total 
charge  delivered  was  found  by  using  the  SEM  efficiency 
factor  of  0.06.  This  represents  the  largest  uncertainty  in 
determining  the  total  charge  delivered.  Due  to  uncertainties 
in  the  efficiency  factor,  the  total  charge  delivered  can  be 
determined  to  no  better  than  5  percent.  Electron  energy  was 
controlled  by  the  use  of  a  magnetic  deflection  system  whose 
stability  and  accuracy  is  within  0.5%.  This  is  considered  to 
be  of  sufficient  stability  and  accuracy  that  any  fluctua- 
tions would  not  be  noticable  within  the  limits  of  this  set 
of  experiments. 

The  off-axis  dose  was  measuned  using  Thermal  Luminescent 
Dosimeters  (TLD's).  TID's  are  a  class  of  radiation  detectors 
that  utilize  an  excited  electron  state  to  register  the  depo- 
sition of  energy.  When  a  TLD  is  heated  the  excited  electron 
returns  to  its  ground  state  and  emits  a  photon.  These 
photons  are  a  direct  measure  cf  the  amount  of  energy  ana 
hence  the  dose  deposited  in  the  material. 
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Finding  a  suitable  container  for  the  liquid  nitrogen 
posed  the  first  major  experimental  problem.  The  TLD's  could 
not  be  subjected  to  cryogenic  temperatures  due  to  possible 
mechanical  failure.  There  was  also  concern  that  whatever  was 
used  to  insulate  the  TLD's  might  affect  the  experiment  in  a 
undesirable   fashion.  Due    to    cost    and      handling   considera- 

tions glass  dewers  were  ruled  out  early  in  the  design  phase. 
Several  methods  for  fabricating  a  container  out  of  commer- 
cially available,  closed  ceil  insulation  were  tried  without 
success.  Finally  a  1  inch  thick,  closed  cell  foam,  sheet 
material  was  found  that  could  be  handled  and  glued  like 
wood.  This  material  was  made  into  boxes  with  interior  cross 
sections  of  20  X  20  cm  and  interior  lengths  of  26,  52,  78 
and    104   cm. 

The  problem  of  where  to  position  the  TLD's  was  solved  by 
finding  the  beam  axis  using  a  phosphor  screen  on  the  TLD 
side  of  the  container  prior  to  filling  with  liguid  Nitrogen. 
This  allowed  positioning  of  the  TLD's  to  an  accuracy  of  0.5 
cm    relative    to    the    beam    axis.  The    TLD's    were    wrapped   in    a 

single  layer  of  0.3  mil  Aluminium  and.  attached  to  the  back 
of  the  container  with  two  sided  tape.  On  several  runs  two 
lines  of  TLD's  were  emplaced,  with  one  set  horizontal  to  the 
ground  and  the  other  vertical;  this  placement  was  to  deter- 
mine the  extent  of  any  beam  asymetries.  Table  IV  list  the 
runs    that    were    made. 
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The  TLD's  were  of  the  Calcium  Fluoride  type  and  had  been 
calibrated  for  use  below  30  MeV.  Extending  the  energy 
observed  to  100  Mev  may  result  in  low  readings;  in  the  event 
there  were  gamma  rays  of  high  enough  energy  that  the  range 
of  the  gamma  ray  produced  particle  is  not  negligible 
compared  to  the  mean  free  path  of  the  primary  particle 
[Ref.  8].  Clearly  this  effect  will  lessen  as  the  length  of 
liquid  nitrogen  is  increased.  Cnf ortunatly,  any  attempt  to 
correct  for  this  under  detection  requires  a  detailed  knowl- 
edge of  the  gamma  spectrum,  which  is  not  available.  Since 
the  shortest  measurement  is  at  a  length  greater  than  0.5 
radiation  lengths  there  will  be  very  few  gamma  photons  with 
energy  greater  than  30  MeV  and  any  effect  of  under  detection 
due  to  high  energy  gamma  photons  should  be  negligible 
[Ref.    1]. 

The  TLD's  consist  of  a  matrix  of  teflon  holding  a  low 
concentration  of  Calcium  Fluoride.  Since  the  absorption 
cross  sections  for  Calcium  Fluoride  differ  by  at  most  5 
percent  from  Silicon  from  0.1  to  100  MeV  the  TLD  measured 
dose    may    be    used    as    dcse   Silicon   [Ref.    9]. 

B.       RESULTS 

All  TLD's  were  read  at  the  Naval  Surface  Weapon  Center, 
White  Oak,  Maryland.  The  TLD's  were  overnight-expressed  to 
NSFC  where  they  were  read  and  the  results  returned  by  mail. 
Results  of  the  experimental  runs  were  reported  in  Eads, 
using  the  charge  delivered  during  the  run,  these  values  were 
converted  to  Rads  per  Coulomb.  Figure  3.1  through  Figure  3.4 
show  the  result  of  the  experimental  runs  at  each  length  of 
liquid  Nitrogen.  Appendix  B  contains  the  complete  data  from 
the    experimental    runs. 

The  vertical  error  bars  in  Figures  3.1  through  3.4  are 
based    on    the   5    percent   non-randcm    uncertainty   in    determining 


the  total  charge  delivered  by  the  Linac  and  the  5  percent 
non-random  uncertainty  resulting  from  the  use  of  Calcium 
Fluoride  TLD's  to  measure  dose  Silicon.  This  results  in  a 
total  uncertainty  of  10  percent  in  the  measured  normalized 
dose.  The  horizontal  error  bars  are  cased  on  the  0.5  centi- 
meter uncertainty  in  the  position  of  the  TLD's  relative  to 
the    learn   axis. 
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Figure   3-1         Measured    Norialized    Dose   at   26    cm. 
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Figure  3.3    Measured  Normalized  Dose  at  78  cm. 
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IV.     CONCLUSIONS 

A.  DCSE    RATE    IN    LIQUID    NITROGEN 

Figures  4.1  through  4.4  show  the  normalized  dose  as 
calculated  by  CYLTEAN  and  measured  with  TLB's.  The  error 
bars  have  been  left  off  of  the  TLD  data  points  in  the 
interest  of  clarity.  Although  there  is  agreement,  within  a 
factor  of  2  at  worst,  between  measured  and  the  calculated 
normalized  doses,  in  many  cases  CYLTRAN  appears  to  calculate 
higher  normalized  doses  than  were  measured.  The  difference 
i^^t'ii!  ttL  tut;  ~c  h  o  nor  iiici  J.J.Z  ed  lo-j-o^  j.  j  uOi.  31'^iiziCc.ii*  jij.  .'  ■■■— y 
it.    due    to    £o:;c    undiscovered    sys  tcr.atic    err.;:. 

iieasurecpnts  at  52,  78,  and  104  cm  and  calculations  ar 
104  cm  show  a  higher  normalized  dose  immediately  off-axis 
than  on-axis.  This  is  due  to  the  small  solid  angle  suDtended 
by  the  on-axis  detector;  any  primary  particle  that  undergoes 
an  interaction  will  be  scatered  into  a  different  detector 
and,  since  the  on-axis  detector  represent  the  smallest  'tar- 
get' of  any  detector,  it  will  have  the  lowest  probability  of 
having    a    particle    scattered    intc   it. 

B.  DOSE    RATE    IN    AIR 

In  order  to  extend  the  above  results  to  ST?  air  it  i;: 
necessary  to  take  Into  account  the  change  in  size,  and 
hence  in  mass,  of  the  detector  as  the  radiation  length 
increases.  The  ratio  of  the  radiaiion  lengths  in  SI?  air  to 
liquid  i-itro^en  is  647  [  Ref .  1].  As  .at  ^tronea  previously 
the  amount  of  energy  in  a  shower  in  some  solid  augle  is  a 
function  of  the  solid  angle  and  the  position  measured  in 
units  of  radiation  length.  If  a  detector  subtends  the  same 
solid  angle  and  is  located  at   the  same  position  as  measured 
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Comparison  of  Experiment  and  Calculations 
at  26  cm. 
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Figure  4.2 
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in  units  of  radiation  length,  then  it  will  absorb  the  same 
amount  of  energy.  However  the  detector  will  be  larger  by  a 
factor  of  the  square  of  the  ratio  of  the  radiation  lengths; 
this  increase  in  size  will  increase  the  mass  and  hence 
decrease  the  dose.  In  the  case  at  hand,  the  dose  absorbed  by 
the  detector  will  decrease  by  647**2.  Table  V  shows  the 
normalized  dose  off-axis  in  STP  air  for  100  MeV  electrons  it 
is  based  on  the  CYTEAN  calculations  in  liquid  Nitrogen,  the 
distance  off-  axis  is  in  terms  of  angle  subtended  by  the 
detector  instead  of  centimeters. 

C.   CCMPAEISON  WITH  PEEYI08JS  PAEERS 

The  only  work  of  a  comparatle  nature  found  in  the  open 
literature  is  [Ref.  10].  This  paper  reports  the  calculated 
off-axis  dose  from  a  50  HeV  electron  beam  in  air  at  ranges 
out  to  3000  meters.  The  calculations  in  [ Eef .  10]  are  based 
on  the  computer  code  ETEAN  which  was  used  to  calculate  the 
electron  and  photon  distribution  emerging  from  a  semi- 
infinite  slab  of  air.  Knowing  the  photon  spectra,  the  dose 
in  a  water  sample  was  calculated.  The  referenced  paper  pres- 
ents the  data  in  Eems  per  electron  with  the  assumption  that 
1  Rem  equals  1  Rad  (H  0). 

Table  VI  compares  selected  data  points  from  the  tables 
in  [Eef.  10]  with  the  corresponding  data  points  from  table  V 
Comparisons  in  this  table  are  made  on  the  assumption  that 
for  that  portion  of  the  dose  due  to  Gamma  photons  1  Rad(Si) 
equals  1  Rad  (HO).  No  account  of  the  difference  in  the 
absorption  cross  sections  for  electrons  needs  to  be  taken 
since  the  referenced  report  assumes  that  no  electrons 
directly  contribute  to  the  dose.  As  can  be  seen  in  table  VI 
the  differences  between  the  normalized  doses  reported  in  the 
referenced  paper  and  the  values  reported  here  are  signifi- 
cant.  At  least  part  of  the   difference  may  be  attributed  to 
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tue  difference  in  incident  energy;  additionally  the  assump- 
tion made  in  [Ref-  10]  that  the  primary,  knock-on  and  the 
photon  produced  (secondary)  electrons  do  not  directly 
contribute  to  the  dose  will  cause  an  under-calculaticn  of 
the  dose.  The  detailed  analysis  necessary  to  determine  the 
exact  source  of  the  differences  is  beyond  the  scope  of  this 
thesis. 


40 


TABLE    V 
Off-Axis    Normalized   Doses   in    STP   Air 


Angular  length    of    STP    air 

displacement     (normalized    dose   in    Hads    per   coulomb    (Si) ) 
(degrees) 

168. 2    m 


0.0-2.2 
2.2-4. 4 
4.4-6.6 
6.6-8.7 
8.7-10.9 
10-  9-13.0 
13-0-15.  1 
15. 1-17. 1 
17. 1-19.1 
19.  1-21.0 

0.0-1.1 
1.1-2.2 
2.2-3.3 
3.3-4.4 
4.4-5.5 
5.5-6.6 
6.6-7.7 
7.7-8.7 
8.7-9.8 
9.8-10.9 


1.856X10* 

1.283X10* 

7.418X103 

3.347X103 

1.250X103 

5.57X102 

2. 11X102 

1.12X102 

6.6X101 

1.8X101 

336. 5    m 
1.617X103 
1.345X103 
1.008X103 
8.07X102 
6.61X102 
5.08X102 
4.06X102 
2.87X102 
2.32X102 
1.75X102 
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T 

Table    V 

Off-Axis 

Normalized    Doses   in    STP    Air    (cont'd.) 

504.  8    m 

0.0-0.7 

5-20X102 

0.7-1.5 

3.77X102 

1.5-2.2 

2.86X102 

2.2-2.9 

2.  12X102 

2.9-3.7 

1.85Z102 

3. 7-4. a 

1. 16X102 

4.4-5.1 

1.09X102 

5.1-5.8 

8.9X101 

5.8-6.6 

6.6X101 

6.6-7.3 

5.0X101 

673.  1    m 

0.0-0.5 

1.39X102 

0.5-1. 1 

1. 71X102 

1.1-1.6 

1.39X102 

1.6-2.2 

1.09X102 

2.2-2.7 

8.7X101 

2.7-3.3 

6.  1X101 

3.3-3.8 

5.7X101 

3.8-4.4 

4.5X101 

4.4-4.9 

3.8X101 

4.9-5.5 

2.9X101 
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TABLE    VI 
Comparison   of   Reported   Results 


Distance 
in   air 
(meters) 

330 

500 
670 


All  values   for    distance    and   normalized    dose    were 
rounded   off    to    the    nearest    decade. 


Angular 
displaceme 
(degrees) 

nt 

CYLIRAN 
Eads/Coul 

Reference 
Rads/Coul 

5.5-10 

220 

20 

0.0-0.5 

520 

60 

0.0-0.5 

140 

30 
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APPENDIX    A 
DESCRIPTION    OF    CY11RAN    AND    OUTPUT 

The  computer  code  CYLTHAN  which  was  used  for  the  calcu- 
lations of  the  off-axis  dose  rates  in  this  thesis  is  a 
cylindrical-geometry  aultimaterial  electron/photon  transport 
code.  The  code  was  developed  at  Sandia  National  Laboratories 
by  J. A.  Halbleib  and  W.H.  Vandevender  [Ref.  11].  The  code 
is  written  in  Fortran  IV  and  was  run  on  a  CDC-7600  at  Los 
Alamos   National    Laboratories. 

CYLTRAN  describes  the  generation  and  the  transport  of 
the  electron/photon  shower  from  the  incident  energy  down  to 
a  cut-off  energy  of  1.0  and  10.0  Kev  for  electrons  and 
photons,  respectively.  The  user  may  specify  a  higher 
cut-off  energy  for  either  or  both  the  electrons  and  photons. 
The  materials  (  up  to  ten  )  that  the  shower  propagates 
through  must  be  invariant  with  respect  to  rotations  about 
some  axis.  The  incident  particles  may  be  either  electrons  or 
photons,  they  may  be  mono-energetic  or  have  a  source  spec- 
trum specified.  The  incident  particles  may  also  be  either; 
mono-directional  (no  divergence) ,  have  a  divergence  angle 
specified  or  be  isotropic.  The  output  of  cyltran  contains 
the    following    main    categories    of   data: 

a)  Charge    and  energy    distribution    profiles. 

b)  Integral  energy  and  number  escape  coefficients 
for    both   electrons   and   photons. 

c)  Escape  coefficients  for  electrons  and  photons 
that  are  differential  in  energy,  angle  and  in  both  energy 
and    angle. 


44 


CY1TRAN  combines  condensed  history  electron  Monte  Carlo 
calculations  with  conventional  single  scattering  Konte  Carlo 
methods-  [Ref.  5]  [fief.  11]-  The  code  follows  electron 
transport  including  energy-loss  straggling,  elastic  scat- 
tering and  the  production  of  knock-on  electrons,  bremsstrah- 
lung,  characteristic  X-rays  and  annihilation  radiation. 
Photon  transport  includes  the  coverage  of  photoelectric, 
Compton  and  pair  production  interactions  and  includes  secon- 
dary   particle    production. 

The  cross  sections  for  electron  interactions  are 
obtained  from  the  ETRAN  Monte  Carlo  system  [Ref.  12]  while 
the  photon  cross  sections  are  calculated  from  analytical 
fits  of  Biggs  and  Lighthill  [fief.  13].  Although  the  problem 
geometry  is  limited  to  two  dinensions  all  particle  geome- 
tries  are    generated    and    followed   in    three    dimensions. 

The  geometry  used  to  generate  the  data  used  here  is 
shown  in  figure  2.1  .  The  problem  geometry  is  rather 
simple,  consisting  of  a  cylinder  of  liquid  nitrogen  ten 
centimeters  in  radius  and  varying  in  length  for  the 
different  runs.  The  Silicon  end  is  set  at  0.2  centimeters  in 
thickness  and  is  divided  into  ten  concentric  rings  each  of 
which  is  one  centimeter  in  width.  Beam  divergence  was  left 
at   zero. 

Only  three  input  files  were  required  to  run  each 
problem  [fief.  11]  table  VII  shows  the  input  parameters, 
their    format      and   their    meaning.  The    problem      inputs    were 

varied  as  to  the  number  of  energy  bins  and  angular  bins  in 
order  to  minimize  the  statistical  uncertainties.  Table  VIII 
is  a  listing  of  the  input  tables  used  in  calculating  the 
off-axis  dose  at  10  4  cm.  The  values  shown  for  the  energy 
bins  and  the  energy  cutoff  values  were  used  for  all  calcula- 
tions. Final  energy  deposition  outputs  for  the  four  problem 
geometries  are  shown  in  table  IX  ,  table  X  ,  table  XI  and 
table    XII 
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APPENDIX  B 
EXPERIMENTAL  RESULTS 

The  tables  in  this  appendix  contain  the  results  of  the 
six  experimental  runs.  Uncertainty  in  charge  delivered  is  5 
percent,  uncertainty  in  the  normalized  dose  is  10  percent 
and  uncertainty  in  off-axis  position  is  0.5  cm. 
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TABLE    XIII 

Ex peri a en tal 

Run 

I 

Length   of   Lie 

Luid 

Nitrogen:       26 

cm. 

Charge   Delivered 

to  SEM: 

1.0    E- 

-06 

COlll. 

Total   Charge 

Delivered: 

1.5    E- 

-05 

coal. 

TLD    POSITION 

DOSE    HfiASIUSfi 

NORMALIZED     DOSE 

(off-axis) 

(Rads) 

(Rads    per    Coul) 

0. 0   cm. 

8.43 

E+04 

5.05 

E+  09 

1.0 

7.00 

E+04 

4.20 

E+09 

2.0 

4.96 

E+04 

2.97 

E+09 

3.0 

3.26 

E+04 

1.95 

E+09 

4.0 

1.75 

E+04 

1.00 

E  +  09 

5.0 

8.  12 

E+03 

4.87 

E+08 

6.0 

5.  13 

E+03 

3.07 

E+08 

7.0 

2.  54 

E+03 

1.52 

E+08 

8.0 

1.32 

E+03 

7.92 

E+07 

9.0 

0.70 

E+03 

4.74 

E+07 

10.0 

... 

0.51 

E+03 

3.06 

E+07 
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.  . 

TABLE    XIV 

Experimental   Run 

II 

Length   of  Liquid 

Nitrogen:       52        cm. 

Charge   Delivered 

to   SEfl: 

4. 4    E-06 

coul. 

Total   Charge    Delivered:    7.33    E-05 

coul. 

TLD    POSITION 

DOSE    MEASUEED 

NORMALIZED    DOSE 

(off-axis) 

(Rads) 

(Rads    per    Coul) 

0.5   cm. 

2.58 

E+04 

3.51 

E+08 

1.5 

4.  41 

E+04 

2.76 

E+08 

2.5 

1.76 

E+0  4 

2.40 

E+08 

3.5 

1.24 

E+04 

1  .69 

E+08 

4.5 

1.03 

E+04 

1.40 

E  +  08 

5.5 

7.9 

E+03 

1  .07 

E+08 

6.5 

6.39 

E+03 

8.71 

E+07 

7.5 

5.29 

E+03 

7.21 

E  +  07 

8.5 

3.83 

E+03 

5.22 

E+07 

9.5 

2.85 

E+03 

3.88 

E+07 
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TABLE    XV 
Experimental   Run   III 

Length    of   Liquid    Nitrogen:       26         cm. 
Charge   Delivered    tc   SEM:     1.00    E-06        coul, 
Total   Charge    Delivered:       1.67    E-05      coul. 


TIP    POSITION 

DOSE    ME AS  DEED 

NORMALIZED    DOSE 

(off- 

axis) 

(Rads) 

(Rads   per    Coul) 

0.0 

cm. 

8.62 

E+04 

5.17 

E+09 

0.5 

8.76 

E+04 

5.25 

E+09 

1.5 

6.83 

E+0<1 

4.09 

E+09 

2.5 

4.  14 

E+04 

2.48 

E  +  09 

3.5 

2.48 

E+04 

1  .47 

E+09 

4.5 

1.38 

E+04 

8.27 

E  +  08 

5.5 

8.  10 

E+03 

4.85 

E+08 

6.5 

4.40 

E+03 

2.64 

E+08 

7.5 

2.  52 

E+03 

1.51 

E+08 

8.5 

1.62 

E+03 

9.69 

E  +  07 

9.5 

1.04 

E+03 

6.24 

E+07 

10.0 

0.52 

E+03 

3.13 

E+07 
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TABLE    2VI 
Experimental   Run   IV 

Length   of   Liquid    Nitrogen:       52        cm. 
Charge   Delivered    to    SEM:     8. 8   E-06        coul. 
Total   Charge    Delivered:     1.35   E-04        coul. 

TIP    POSITION  IOSE    MEASUjvED  NORMALIZED    DOSE 

(off-axis)  (Rads)  (Rads    per    Coul) 


5. 17  E+08 
5.53    E+08 

4.82  E+08 
3.42  E+08 
2.77  E+08 
2.28  E+08 

1.83  S+08 
1.46  E+08 
1.20  E+08 
8.07   E+07 


(horizontal) 

0.0 

cm. 

6.99 

E+04 

0.5 

7.47 

E+04 

1.5 

6.51 

E+0  4 

2.5 

4.62 

E+0  4 

3.5 

3.74 

E+04 

4.5 

3.08 

E+04 

5.5 

2.47 

E+04 

6.5 

1.97 

E+04 

7.5 

1.62 

E+04 

8.5 

1.09 

E+04 

9.5 

failed 
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Table  XVI 

Exp 

erimental 

Run  IV 

(cont'd.) 

TIP  POSITION 

DOSE  MEASURED 

NORMALIZED  DOSE 

(Off- 

axis) 

(Rads) 

(Rads  per  Coul) 

(ver 

tical) 

0.5 

5.  94 

E+04 

4.4.0 

E+08 

1.5 

4.90 

E  +  04 

3.62 

E+08 

2.5 

failed 

3.5 

3.79 

E  +  04 

2.80 

E+08 

4.5 

2.94 

E+04 

2.  17 

E+08 

5.5 

2.40 

E+0  4 

1.77 

E+08 

6.5 

1.81 

E+04 

1.34 

E+08 

7.5 

1.57 

E+04 

1.16 

E+07 

8.5 

1.17 

E+0  4 

8.67 

E+07 

9.5 

7.85 

E  +  03 

5.81 

E+07 

10.5 

6.93 

E+03 

5.  13 

E+07 

, ._.  . 

i 

65 


r  " 

T 

TABLE  X7II 

Experimental 

Run 

7 

Length  of  Liquid 

Nitrogen:   78 

cm 

• 

Charge  Delivered 

to  SEM: 

2.8  E- 

05 

coul. 

Total  Charge  Delivered:  4.75  E- 

■04 

coul. 

TID  POSITION 

DOSE  MEASURED 

NORMALIZED  DOSE 

(off-axis) 

(Rads) 

(Rads  per  Coul) 

(horizontal) 

0.0  cm. 

8.  63 

E+04 

1.81 

E+06 

0.5 

8.97 

E+04 

1.88 

E+08 

1.5 

7.  19 

E+04 

1.51 

E+08 

2.5 

5.06 

E+04 

1.06 

E  +  08 

3.5 

3.  76 

E+04 

7.92 

E+07 

4.5 

3.  15 

E+04 

6.63 

E+07 

5.5 

2.37 

E+04 

4.99 

E  +  07 

6.5 

1.83 

E+0  4 

3.85 

E  +  07 

7.5 

1.50 

E+04 

3.  16 

E+07 

8.5 

1.22 

E+0  4 

2.57 

E+07 

9.5 

9.31 

E+03 

1.96 

E+07 

10.0 

9.63 

E+03 

2.03 

E+07 

i 

j 
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■ 

table  XVII 

Exp 

erimental  Run  V 

(cont'd.) 

TIP  POSITION 

DOSE  MEASOEED 

NORMALIZED  DOSE 

(off-axis 

) 

(Rads) 

(Rads  per  Coul) 

{vertica 

1) 

0.5 

7.33 

E+OU 

1.54 

E+08 

1.5 

6.95 

E+04 

1.46 

E+08 

2.5 

4.27 

E+0  4 

8.99 

E  +  07 

3.5 

3.68 

E+0  4 

7.75 

E+07 

4.5 

2.60 

E+04 

5.47 

E+07 

5.5 

2.00 

E  +  04 

4.21 

E+07 

6.5 

1.53 

E+OU 

3.22 

E+07 

7.5 

1.  17 

E+04 

2.46 

E+07 

8.5 

1.02 

E+04 

2.  15 

E+07 

9.5 

8.03 

E+03 

1,69 

E+07 
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I 

i 

TABLE  XVIII 

Experimental  Ran 

71 

Length  of  Liquid  Nitrogen:   104   cm. 

Charge  Delivered  tc  SEM: 

5.0  E-0  5 

coul. 

Total  Charge  Delivered:  8.33  E-04 

coul. 

TLD  POSITION      LOSE  MEASURED 

NORMALIZED  DOSE 

(off-axis)         (Rads) 

(Rads  per  Coul) 

(horizontal) 

0.0  cm.         5.98 

E+04 

7.18 

E+07 

0.5             6.13 

E+04 

7.36 

E+07 

1.5             5.46 

E+04 

6.55 

E+07 

2.5             4.48 

E+04 

5.38 

E  +  07 

3.5             3.93 

E+04 

4.72 

E+07 

4.5             3.03 

E+04 

3.64 

E+07 

5.5             2.51 

E+04 

3.01 

E+07 

6.5             1.99 

E+04 

2.39 

E  +  07 

7.5             1.63 

E+04 

1.96 

E+07 

8.5             1-24 

E+04 

1.49 

E  +  07 

9.5             1.04 

E+04 

1.25 

E+07 

10.0             9.55 

E+03 

1.  15 

E+07 

j 
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Table  XUIII 

Exp 

erimental 

Run  VI 

(cont'd.) 

TID  POSITI 

ON 

DOSE  MEASUKED 

NORMALIZED  DOSS 

(off- 

axis) 

(Bads) 

(Rads  per  Coul) 

(ver 

ti 

cal) 

0.5 

6.  13 

E+04 

7.36 

E  +  07 

1.5 

5.42 

E+0  4 

6.51 

E+07 

2.5 

3.87 

E+04 

4.65 

E+07 

3.5 

3.21 

E+04 

3.85 

E+07 

a.  5 

2.62 

E+04 

3.  14 

E+07 

5.5 

2.  14 

E+04 

2.57 

E+07 

6.5 

1.66 

E+04 

1.99 

E  +  07 

7.5 

1.34 

E+04 

1.61 

E+07 

8.5 

1.08 

E+04 

1.30 

E+07 

9.5 

8.70 

E+03 

1  .04 

E+07 
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c-  1        Comparison  of  calcu- 
lations and  measure- 
ments of  the  off-axis 
radiation  dose  (Si)  in 
liquid  nitrogen  as  a 
function  of  radiation 
length. 


